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INTRODUCTION 


We  have  employed  BENCHMARK-76,  a version  of  the  AIRCHEM  code  with 
a 64  species  set  and  a nominal  496  reaction  set,  to  obtain  species 
densities  as  a function  of  time  at  an  altitude  of  60  km  for  the  foll^Y^'^S 
conditions^  For  the  prompt  ionization,  N^  was  assigned  the  values  10  , 

10^  or  10°  ion-pairs  cm  . The  delayed  ionization  is  given  by 

Q{t)  = Qq  (1  + (1) 

6 8 

where  t is  the  time  in  seconds,  and  assigned  the  values  10  , 10  or 
10^®  ion-pairs  cm”^  s"^,  subject  to  the  constraint  |q^(  < 1n^|.  The 

relative  distribution  of  the  species:  N_*,  0*,  N*,  0^,  N(^S),  N(^D) , 

3 11  1 ^ ^ 

0(  P) , 0(  D)  02(a  A)  and  02(b  Z)  produced  by  ionization  under  disturbed 

conditions  is  taken  to  be:  0.64,  0.16,  0.14,  0.06,  0.45,  0.61,  1.28, 

0.10,  0.24  and  0.05  particles  per  ion  pair^  Respectively.  This  dis- 
tribution closely  follows  that  of  Gilmore,  ’ except  that  N2(A^Z)  is 
not  carried,  [02(^A)]  is  taken  as  0.24  and  [0]  is  1.28.  Cases  for  daytime 
and  nighttime  conditions  are  reported;  the  calculations  correspond  to 
about  a three  hour  interval  centered  around  noon  and  ^idnight.  Approp- 
riate neutral  densities  have  been  outlined  elsewhere.  The  program  coding 
and  fixed  input  parameters  have  also  been  discussed  previously. 

It  is  the  purpose  of  this  report  to:  (1)  display  selected  computed 
results  as  a function  of  time  and  ionization  conditions,  (2)  show  the 
results  of  comparison  studies  of  a very  limited  nature,  and  (3)  report 
the  computed  variations  of  the  equivalent  rate  coefficients  (ERC's)  as  a 
function  of  time  and  ionization  conditions.  (A  copy  of  these  BENCHMARK- 
76  results  may  be  obtained  by  writing  the  authors.) 


^E.  L.  Lortie,  M.  D.  Kregel  and  F.  E.  Niles,  "AIRCHEM:  A Computational 
Technique  for  Modeling  the  Chemistry  of  the  Atmosphere,"  BRL  Report 
No.  1913,  August  1976.  (AD#  A030157) 

2 

F.  Gilmore  as  quoted  by  B.  F.  Myers  and  M.  R.  Schoonover,  "Electron 
Energy  Degradation  in  the  Atmosphere:  Consequent  Species  and  Energy 
Densities,  Electron  Flux,  and  Radiation  Spectra,"  DNA  3513T,  3 Jan  75, 
Table  6. 

^J.  M.  Heimerl  and  F.  E.  Niles,  "BENCHMARK- 76:  Model  Computations 
for  Disturbed  Atmospheric  Conditions  I.  Input  Parameters,"  BRL 
Report  No.  2022,  October  1977.  (AD  #A050355) 

4 

F.  E.  Niles  and  J.  M.  Heimerl,  "Selected  Neutral  Species  Profiles 
0-100  km,"  BRL  Memo  Report  2767,  July  1977.  (AD#  A042620) 
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RESULTS 


A.  Selected  Densities 

The  daytime  and  nighttime  electron  density  computed  as  a function 
of  time  after  burst  is  shown  for  the  six  ionization  conditions  in  Figures 
1 and  2,  respectively. 

We  might  expect  that  the  larger  the  disturbance,  i.e.,  the  larger 
the  value  of  N for  a given  Q , the  longer  it  would  take  for  the  electron 
density  to  return  to  its  quiescent  value.  _In  Figure  1 curves  4 and  6 
(also  2 and  3)  at  the  latest  times  (t  > lO'^s)  show  that  the  computed  value 
of  the  electron  density,  [e] , for  curve  4 (2)  falls  below  that  for  curve 
6 (3) . To  examine  why  our  expectations  are  not  borne  out  by  these  compu- 
tations, let  us  focus  attention  on  curves  4 and  6 of  Figure  1. 

For  both  of  these  curves  the  main  charge  flow  is  found  to  be 
Q ->•  e -»■  O2"  COj".  The  major  reactions  are 

e + 0^  + O2  ->  02'  + O2  R125* 

O2"  + 03-^  O3"  + ^2  R246 

and  Oj  CO2  CO^  0.,.  R269 

A secondary  path  exists;  to  wit, 

0^'  + 0,  + M -V  0.'  + M R285 

2 2 4 

0/  + CO^  CO."  M R273 

4 2 4 

and  CO.  + 0 -*•  C0_  + 0, . R255 

4 3 2 

In  any  event  the  charge  at  this  point  in  the  sequence  resides  upon  the 
COj"  ion.  The  major  loss  process  for  this  ion  is  found  to  be  the  reaction 

COj"  + 0 ->■  02"  + CO2.  R252 

which  simply  cycles  charge  back  to  the  O2  ion.  The  major  loss  of  process 
for  the  O2”  ion  in  both  cases  is  the  reaction 

O2'  + 02C^A)  -»■  2 O2  + e,  R137 

whose  product  is  an  electron  and  the  entire  cycle  e j ©2"  CO^" 
begins  again. 

For  curve  4 there  exists  an  important  secondary  loss  process  for 
the  COj"  ion;  namely,  the  reaction, 

COj"  + NO  -*■  NO2"  + CO2.  R250 

* — 

Read  as  reaction  125,  see  Appendix  for  complete  reaction  listing. 
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Figure  1.  Logarithm  of  the  computed  daytime  electron  density  at  60  km 
as  a function  of  the  logarithm  of  time  for  six  different 
excitation  conditions.  Q is  the  delayed  ionization  parameter  ^ 
in  ion-pairs  cm'^s"  and  the  prompt  ionization  parameter  in  cm""^. 
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Figure  2.  Same  as  Figure  1,  except  for  nighttime  conditions. 


L 


10 


This  reaction  allows  about  10%  of  the  charge  to  flow  to  the  NO2  ion 
which  reacts  rapidly  to  form  the  NO^  ion  by  the  reaction 

N02“  + 03-^  NO^’  02‘  ^ 

By  contrast,  reaction  250  for  the  conditions  of  curve  6 allows  only 
about  0.1%  of  the  charge  to  flow  to  the  NO^  ion.  It  is  the  relative 
importance  of  the  CO^'  ion  loss  via  reaction  250  that  accounts  for  the 
calculated  difference  in  the  late  time  electron  densities.  For  case 
4 the  charge  that  flows  to  the  NO^'  ion  tends  to  remain  bound  to  it. 

(The  electron  affinity  of  the  ion  (EA  = 3.9  eV)  is  greater  than 

that  of  the  O2  ion  (0.440  eV)  or  the  C0_  ion  (EA  = 2.69  eV)  .) 

This  process  tends  to  deplete  the  number  density  of  free  electrons. 

For  curve  6,  the  path  forming  the  NO^  ion  is  about  a hundred  times  less 
effective  than  for  curve  4 and  the  number  density  of  free  electrons  tends 
to  be  maintained  through  the  recycling  process. 

To  find  why  reaction  250  is  so  pivotal  a reaction  the  values  of  [NO] 
are  examined  at  late  times.  It  is  found  that  the  value  of  [NO]  for  curve 
4 is  about  one  hundred  times  larger  than  for  curve  6.  The  reason  for 
this  lies  in  the  formation  of  the  NO  molecule  at  the  earliest  times  by 
the  reaction 

N(^D)  + O2  NO  + 0.  R86 

2 

The  value  of  [N(  D)],  hence  the  value  of  [NO]  through  reaction  86,  is 
determined  in  large  part  by  the  value  of  the  prompt  ionization  parameter, 
N . To  summarize  we  find  that  for  the  same  value  of  Q,  the  greater  the 
value  of  N^,  the  larger  the  value  of  [NO],  the  larger  the  value  of 
[NOj  ] , the  smaller  the  value  of  [e] . Similar  arguments  apply  for 
curves  2 and  3 of  Figure  1. 

Figure  2 explicitly  shows  a similar  phenomenon  at  the  latest  times 
for  curves  2 and  3.  Curves  4 and  6 also  cross  as  in  Figure  1 but  their 
values  are  so  similar  to  case  5 that  they  are  indistinguishable  on  this 
scale.  Arguments  similar  to  the  daytime  case  apply  here,  except  that 
the  details  of  the  reactions  change.  As  before,  the  main  charge  flow 
is  given  by  Q ->■  e O2  -*•  CO^”  by  reactions  125,  246  and  269.  But  here 
the  major  loss  for  the  CO^  ion  at  late  times  at  night  is  given  by 

COj"  + NO2  NOj"  + CO2.  R 251 

^See  for  example  E.  E.  Ferguson,  D.  B.  Dunkin  and  F.  C.  Fehsenfeld, 
"Reactions  of  NO^”  and  N0_~  with  HCl  and  HBr,"  J.  Chem.  Phys. 

1459-1463,  1972. 

^R.  J.  Celotta,  R.  A.  Bennett,  J.  L.  Hall,  M.  W.  Siegel  and  J.  Levine, 
"Molecular  Photodetachment  Spectrometry.  II.  The  Electron  Affinity 

^of  O2  and  the  Structure  of  02',..  phys.  Rev.  631-642,  1972. 

^S.  P.  Hong,  S.  B.  Woo  and  E.  M.  Helmy,  "Photodetachment  of  Thermally 
Relaxed  CO^',"  Phys.  Rev.  A15,  1563-1569,  1977. 
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The  NOj'  ion  formation  is  aided  among  other  reactions  by  the  secondary 
reactions 

02'  + NO2  -<■  N02‘  + O2  R244 

followed  by 

N02’  +03-^  NO^’  + 0^.  R260 

There  are  two  differences  for  nighttime  conditions.  1)  The  NO-  molecule 
rather  than  the  NO  molecule  determines  the  charge  flow  to  the  NO,”  ion 
and  2)  the  recycling  of  charge  is  not  important  at  late  times.  '*This 
last  difference  comes  about  because  reaction  252  is  not  an  important  loss 
process  for  CO^  at  night.  The  fundamental  reason  for  this  lies  in  the 
fact  that  the  value  of  [0]  falls  markedly  at  late  times  at  night  relative 
to  the  daytime  case.  Atomic  oxygen  is  converted  to  ozone  by  the  reaction 

0 + O2  + M ->■  0^  + M;  R117 

but  the  dominant  daytime  loss  process  for  ozone,  the  reaction 

O5  + hv  ->■  0(^D)  + 02(^A)  , R47 

which  returns  atomic  oxygen  through  the  reactions 

OC^O)  + M ^ 0 + M R97 

and  0(^D)  + 0^  ->■  0 + , R99 

is  unavailable  at  night.  A similar  photon  reaction  explains  the  NO/NO2 
difference.  To  wit,  at  late  times  the  primary  source  of  the  NO2  molecule 
is  the  reaction 


NO  + Oj  -»■  NO  + O2.  R88 

In  daytime  the  primary  loss  process  for  the  NO2  molecule  is  the  reaction 

NO2  + hv  ->  NO  + 0 R40 

with  aid  from  the  secondary  reaction 

NO2  + 0 NO  + O2.  R92 

Again,  at  night  when  no  photons  are  available  reaction  40  is  not  signifi- 
cant. Reaction  92  does  not  compensate  for  the  loss  because  as  we  have 
alluded  to  above  the  atomic  oxygen  density  also  falls  at  late  times  at 
night  with  the  result  that  the  [NO2]  builds. 
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The  dominant  electron  depopulating  process  at  late  times  for  both 
day  and  night  conditions  (Fig’jires  1 and  2)  is  three  body  attachment. 

At  night  the  dominant  electron  production  mechanism  at  late  times  is 
the  delayed  ionization  source  term  whose  value  is  given  by  Q (see  equa- 
tion 1) . Thus  at  the  later  times  the  electron  densities  tend  to  track 
Q.  For  daytime  conditions  the  recycling  of  charge,  as  discussed  above, 
may  permit  the  rates  of  other  reactions  to  compete  with  the  source  rate 
with  the  result  that  the  daytime  electron  densities  need  not,  in  general, 
track  Q. 

Figures  3 and  4 show  the  sum  of  the  negative  ions  (exclusive  of 
electrons)  computed  as  a function  of  time.  Conditions  correspond  to 
those  in  Figures  1 and  2,  respectively.  Figures  5 and  6 show  the  corre- 
sponding plots  for  the  sum  of  the  positive  ions. 

-2 

For  times  less  than  10  seconds,  the  electron  and  positive  ion 
curves  figures  1,  2,  5 and  6)  tend  to  group  according  to  the  prompt 
ionization  parameter,  N . For  times  greater  than  10*  seconds  they  tend 
to  group  according  to  the  delayed  ionization  parameter,  Q^. 

Some  of  the  curves  cease  abruptly  in  these  (and  other)  figures 
because  computations  are  automatically  halted  whenever  the  electron 
density  falls  below  10  cm  , the  assumed  quiescent  background  electron 
density. 


B.  Comparison 

Figures  7 and  8 show  the  electron  density  computed  as  a function  of 
time  at  60  km  altitude  fgr  both  the  present  results  (BRL-76)  and  the 
DCHEM  results  of  Scheibe  for  daytime  and  nighttime,  respectively.  In 
spite  of  the  fact  that  different  sets  of  chemical  reactions  and  inte- 
gration schemes  were  employed,  the  agreement  is  very  good  to  ex^gllent. 
Both  day  and  night  comparisons  have  alsogbeen  made  and  jepo^ted  for 
the  fix|^  ionization  conditions  Q = 10°  ion-pairs-cm  -s  and 
N = 10  cm  . These  comparisons^indicate  a basic  overall  agreement 
between  Scheibe' s results  and  our  own  for  the  computation  of  the  electron 
density  at  60  km  as  a function  of  Day/Night  conditions,  of  excitation 
and  of  time. 


g 

W.  S.  Knapp,  "A  Simplified  D-Region  Chemistry  Model  for  Nuclear 
Environments,"  DNA  2850T,  April  72;  also  M.  Scheibe,  private  communi- 
cation, 1976. 

9 

For  example,  the  BENCHMARK-76  reaction  set  includes  neutral  odd-hydrogen 
reactions;  Scheibe' s set  does  not. 

M.  Heimerl  and  F.  E.  Niles,  "BENCHMARK-76:  Model  Computations 
for  Disturbed  Atmospheric  Conditions  II.  Results  for  the  Stratosphere 
and  Mesosphere,"  BRL  Technical  Report  ARBRL-TR-02050,  March  1978. 
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Figure  3.  Logarithm  of  the  daytime  total  negative  ion  density  60  km  as 

a function  of  the  logarithm  of  time  for  six  different  ionization 
conditions. 
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Figure  4.  Same  as  Figure  3,  except  for  nighttime  conditions. 
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Figure  5.  Logarithm  of  the  daytime  total  positive  ion  density  at  60  km  as 

a function  of  the  logarithm  of  time  for  six  ionization  conditions. 
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Figure  7.  For  two  excitation  conditions  at  60  km  the  logarithm  of  the 
daytime  electron  density  as  computed  by  BRL's  BENCHMARK-76 
code  H and  Schiebe's  DAIRCHEM  code  (X)  vs.  the  logarithm 
of  time. 
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Figure  8.  Same  as  Figure  7,  except  for  nighttime  conditions. 
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C.  Equivalent  Rate  Coefficients  (ERC's) 

Once  the  individual  species  concentrations  have  been  computed  as  a 
function  of  time  and  excitation  conditions  a set  of  ERC's  can  be  gen- 
erated. 

The  effective  electron-ion  recombination  coefficient  is  defined  by 

= 5:Rg./[e]E[N^^]  (2) 

where  the  R . are  the  rajes  for  electron- ion  recombination,*  [e]  the 
electron  density  and  [Nj^  ] is  the  density  of  the  positive  ion  species. 
The  sums  are  over  all  positive  ion  species. 

Figures  9 and  10  show  the  temporal  variation  of  a,  for  daytime  and 
nighttime  conditions,  respectively.  In  both  of  these  figures,  ranges 
from  values  characteristic  of  electron-ion  recombination  with  diatomic 
ions  at  early  times  to  values  characteristic  of  recombination  with  clus- 
tered ions  at  late  times.  The  return  to  the  quiescent  clustered  ion 
condition  seems  to  be  dominated  by  the  value  of  Q^,  the  delayed  ioniza- 
tion pagameter.  To  wig,  jn  curves  4,  5 and  6,  corresponding  to 

= 10°  ion-pairs  cm"'*s'  , the  cluster  ions  seem  to  be  reestablished 

after  about  100  seconds.  In  curves  2 and  3,  corresponding  to 

= 10®  ion-pairs  cm~^s'\  this  time  is  greater  than  10^  seconds  during 

the  day  (Figure  9)  and  about  10^  seconds  at  night  (Figure  10) . For 

curve  1 (0^  = 10^^  ion-pairs  cm’^s”^)  this  time  appears  to  be  even 
longer. 

The  effective  electron  attachment  frequency  is  given  by 

A = ZR^/[e]  (3) 

where  the  R.  are  the  rates  for  attachment  processes  and  [e]  is  the  electron 
density.  The  sum  is  over  all  attachment  processes.** 

Figures  11  and  12  show  A vs.  time  for  daytime  and  nighttime  condi- 
tions, respectively.  Different  curves  correspond  to  differing  ionization 
conditions.  The  variations  in  these  curves  will  now  be  traced  to  varia- 
tions in  the  value  of  [0^] . 

From  the  output  we  find  that  ZR^  = Rj24  "*■  ^^25  *^126 

Oj  + e ->  0"  + O2  R124 

©2  + e + ©2  O2  * O2  R125 

^2  ♦ e + -*■  ^ ^2  R126 

V 

Reactions  138-162,  inclusive  (see  Appendix) 

** 

Reactions  122-126,  inclusive  (see  Appendix) 
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Figure  9. 


Logarithm  of  the  computed  daytime  effective  electron-ion  recombi 
nation  coefficient,  at  60  km  as  a function  of  the  logarithm 
of  time  for  six  ionization  conditions. 
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Figure  11.  Logarithm  of  the  computed  daytime  effective  attachment  frequency, 
A,  at  60  km  as  a function  of  the  logarithm  of  time  for  six 
ionization  conditions. 
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log  r (s) 

Figure  12,  Same  as  Figure  11,  except  for  nighttime  conditions. 
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Substituting  into  equation  (3)  we  find  that 

A = C(1  + B)  (4) 

where  C = k^25 [^2^ ^ ^ *^126 

and  B = ^^241033/0.  (6) 

Since  C in  Equation  (5)  is  constant  to  several  decimal  places  and  since 
kj24  in  Equation  (6)  is  held  fixed,  the  changes  in  the  value  of  A in 
Equation  (4)  must  reflect  changes  in  the  ozone  concentration. 


For  the  daytime  case.  Figure  11,  let  us  consider  curve  4 since  it 
has  the  largest  excursion.  Figure  11  shows  three  easily  determined  time 
-3  -1 

domains.  For  10  < t < 10  s A is  observed  to  be  constant.  It  is 

found  that  B of  Equation  (6)  is  0.01  in  this  time  period.  [B  of  Equation 
(6)  is  simply  a measure  of  the  relative  importance  of  ^124*^  ^*^125  ^ *^126^'^ 


Thus  Equations  (4)  and  (5)  show  that  at  early  times  A » C.  The  reason 
that  the  value  of  [0^]  remains  constant  over  this  interval  lies  in  the 
near  balance  between  the  production  and  loss  rates  of  ozone.  The  major 
reactions  populating  ozone  are 


0 + O2  + M ->■  0^  + M R117 

and  0^^  + 0 -►  02^  + 0^.  R204 

Reacting  204  is  found  to  be  an  important  source  of  0^  up  to  and  including 
t = 10  s.  The  major  reaction  depopulating  ozone  over  this  time  interval 
is 


02(^1)  + O5  2 O2  + 0.  R104 

The  overall  effect  of  reactions  117  and  104  is  to  convert  02^^P  ®2’ 

Eventually  [O-C^S)]  is  depleted  and  for  times  greater  than  10  ^s  the  value 
of  [(02(^5;)]  rails  to  such  an  extent  that  the  production  rate  (now  due 
almost  exclusively  to  Reaction  117)  exceeds  the  l^ss  rate  with  the 
result  that  the  ozone  concentration  rises.  By  10  seconds  the  ozone 
concentration  has  risen  to  such  an  extent  that  the  rate  of  0^  loss  by 
the  reaction 

Oj  + hv  ->■  0(^D)  + 02(^A)  R49 

3 

tends  to  balance  the  0_  production  rate.  For  times  greater  than  10  s 
the  value  of  [0-]  tracks  the  value  of  [0]  (in  fact  their  densities  are 
about  equal)  and  since  the  value  of  the  atomic  oxygen  begins  to  fall, 
the  value  of  the  ozone  density  does  also. 


At  night  only  two  time  domains  are  observed  in  Figure  12.  For 
times  less  than  10°s  Jhe  daytime  arguments  apply  for  A constant.  For 
times  greater  than  lO'^  seconds  the  loss  rates  which  are  due  to 


H + Oj  ->■  HO  + O2  R70 

N + Oj  NO  + O2  R85 

NO  + Oj  ->■  NO2  + O2  R88 

0 + Oj  ->■  O2  + O2  R93 

0 + Oj  -*■  02(^A)  + O2  R94 


and  Reaction  (104)  fall  relative  to  the  production  by  Reaction  (117) 
and  again  the  rise  in  the  value  of  [0^]  is  reflected  by  an  increase  in 
the  value  of  A in  Figure  12. ^ The  ozone  density,  [0^],  will  not  rise 
indefinitely.  In  fact  at  lO'^s  the  major  loss  rate;  which  is  due  to 
Reaction  88,  greatly  exceeds  the  production  rate  to  117  thus  we  would 
anticipate  the  value  of  [0_]  (hence  the  value  of  A)  to  fall  as  in  the 
daytime  case  but  at  times  later  than  those  graphed  in  Figure  12. 

The  effective  electron  detachment  frequency  is  defined  by 

D = ZRp/E[N^'],  (7) 

where  the  R^  a^^  the  rates  for  individual  detachment  processes*  and 
[N, ']  is  the  negative  ion  density.  The  sums  are  over  all  detachment 
rates  and  all  negative  ion  densities,  respectively. 

Figures  13  and  14  show  D vs.  time  for  daytime  and  nighttime  condi- 
tions, respectively.  Ionization  conditions  are  listed  nearby  their 
corresponding  curves.  The  most  notable  feature  in  both  Figures  13  and 
14  is  the  apparent  grouping  of  the  curves.  There  are  three  groups, 
the  first  composed  of  curves  1,  2 and  4;  the  second  of  curves  3 and  5 
and  the  third  of  curve  6.  From  the  key  on  these  figures  we  find  that 
they  are  grouped  according  to  the  prompt  ionization  parameter,  N^. 

In  Figure  13  curve  6 attains  the  approximate  limit  of  10°s"^  after 
10  seconds.  This  limiting  value  is  typical  of  daytime  quiescent  con- 
ditions and  is  the  value  the  other  five  daytime  curves  should  approach, 
but  do  not.  We  shall  now  address  (1)  why  curve  6 attains  its  limiting 
value  and  (2)  why  the  others  do  not  attain  a similar  limit.  To  be 
specific  we  shall  concentrate  our  attention  on  curve  4 as  a contrast 
to  curve  6. 

It  is  found  to  a good  approximation  that 


^*^0  * '^137 


*^135 


129’ 


Reactions  1-5  and  11  for  photodetachment;  127-137  inclusive  for  two  body 
collisional  detachment.  See  Appendix  for  reaction  list. 
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Figure  14.  Same  as  Figure  13,  except  for  nighttime  conditions. 
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substituting  this  expression  into  Equation  (7)  we  have 
D«  (k^3^[02(^A)]  + kj35[0])  ([02’]/I:[N'])  + kj29[N2] [0']/Z[N^'] . (8) 

For  case  6 we  see  from  Figure  15  that  for  times  greater  than  10°s 
^[\"]  **  addition  [02~]/[C03"]  and  [0”]/[C03']  are  approxi- 

mately constant  over  this  time  interval  (see  Fig.  15).  Since  the  value 
of  (N2]  and  the  rate  coefficients  are  constant  the  behavior  of  the 

0~(^A)  molecule  and  of  the  0 atom  needs  to  be  examined.  Over  the  interval 
^14  1 

10  < t < 10  s,  the  value  of  [O-C^A)]  monotomically  rises  from  1.5  x 
10  3 10  -3  ^ 

10  cm  to  3.4  X 10  cm  . There  is  no  dramatic  change  in  the  value  of 

[0]  either.  It  first  declines  from  1.4  x lO^^cm”^  to  1.1  x lO^^cm 

10  -3 

then  rises  to  1.3  x 10  cm  . Since  k,__  and  k..-  are  comparable 
-10  3 -10  3 

(2  X 10  cm  /s  vs.  1.5  X 10  cm  /s,  respectively)  the  variations  in 
the  values  of  [02(^A)]  and  [0]  tend  to  offset  each  other.  Thus  for  case 
6 and  for  times  greater  than  10°s  D is  approximately  constant. 


On  the  other  hand,  Figure  16  shows  for  case  4 and  for  times  greater 
than  lO^s  that  E[Nj^  ] « [NO^  ].  In  addition  both  [O2  ]/[N03  ] and 

[0  ]/[N0-"]  can  be  seen  from  the  figure  to  be  declining  functions  of 

^ 11-3 

time.  The  value  of  [0]  monotomically  declines  from  4.4  x 10  cm  to 

1.0  X lO^^cm”^  (10^  < t < lO^s),  while  the  value  of  [02(^6)]  first  rises 

from  4.0  x lO^^cm^  (t  = lO^s)  to  3.4  x lO^^cm  ^ (t  = lO^s) , then  falls 

to  2.7  X lO^^cm”^  (t  = lO^s) . Since  k^^^  *^135  comparable,  the 

declining  value  of  [0]  tends  to  offset  the  rising  value  of  [02(^A)]. 

Thus  D in  Equation  (8)  will  follow  the  ratio  of  the  negative  ion  den- 
sities ^nd  the  net  result  is  that  D for  case  4 and  for  times  greater 
than  10  s is  a decreasing  function  of  time.  Cases  1,  2,  3 and  5 are 
similar  to  case  4. 

One  can  ask  why  the  late  time  dominant  negative  ion  is  NO.”  in 
case  4 (and  in  other  cases)  but  not  in  case  6.  The  answer  to  this  is 
found  to  be  in  the  production  of  the  NO_'  ion  from  either  the  CO."  or 


3^  ions.  The  specific  reactions  are 


^ound 


to  be 


followed  by 


CO.  + NO  -*■  OONO  + CO. 
4 2 


OONO"  + NO  -►  NO  ■ + NO 
2 2 
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Figure  15.  Histories  of  selected  daytime  negative  ion  densities  at  60  km 
for  thegexcitatlon  conditions:  _Q  = 10®  Ion-pairs  cm'^s"'  and 
N = 10  cm  (Case  6).  The  NO^  Ton  density,  dashed  line.  Is 
shown  for  easy  comparison  with  Case  4,  Fig.  16. 
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and 

COj"  + NO  -►  NO2"  + CO2.  R250 

In  either  case  the  NO^  ion  is  formed  from  the  NO2"  ion  by  the 
reaction 

^^2"  °3  ^°3'  °2* 

Reactions  250  and  254  are  minor  loss  processes  for  the  CO_”  and  CO."  ions, 
respectively;  and,  together  with  Reaction  262  depend  critically  on^the 
value  of  [NO] . Figure  17  shows  [NO]  as  a function  of  time  for  all  six 

daytime  cases.  In  case  6 the  value  of  [NO]  rises  to  about  a factor  of 

three  higher  than  the  quiescent  case  (dashed  line)  but  in  case  4 the 
value  of  [NO]  rises  more  than  a factor  of  100  above  the  quiescent  case. 
These  order  of  magnitude  changes  in  the  value  of  [NO]  relative  to  the 
quiescent  case  enable  the  above  reactions  to  enhance  the  NO-  ion 
production  relative  to  case  6. 

Figure  17  also  shows  the  correlation  between  the  value  of  [NO]  and 
the  primary  ionization  parameter  N^.  This  correlation  can  be  explained 
by  considering  the  reactions  producing  the  NO  molecule  at  the  earliest 
times.  They  are  in  order  of  importance 

N(^D)  + O2  ->■  NO  + 0,  R86 

N + O2  -►  NO  + 0,  R84 

and  N + Oj  -*■  NO  + O2.  R85 

The  early  time  values  of  the  atomic  nitrogen  densities  correlate  with 
the  prompt  ionization  parameter,  N . (See  Introduction  for  the  relative 
distribution  of  species  produced  by  ionization  under  disturbed  conditions.) 

At  the  later  times  in  Figure  17  the  value  of  [NO]  for  a given  set 
of  excitation  conditions  is  approximately  constant  and  the  model  results 
predict  that  a steady  state  condition  exists  between  the  NO  and  NO- 
molecules.  At  the  later  times  the  NO2  molecule  is  produced  from  the 
NO  molecule  primarily  by  the  reaction'^ 

NO  + Oj  -►  NO2  + O2  R88 

and  secondarily  by  the  reaction 

NO  + 0 + M -*•  NO2  + M,  R115 

while  the  NO  molecule  is  produced  from  the  NO2  molecule  primarily  by 
the  reaction 

NO2  + hv  ->•  NO  + 0 R40 
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and  secondarily  by  the  reaction 


NO2  + 0 -*■  NO  + O2.  R92 

These  four  reactions  are  the  major  production  and  loss  reactions  for 
both  molecules.  Other  computations  have  exhibited  this  coupling  between 
the  NO  and  NO2  molecules.^ 

Could  the  constant  values  of  [NO]  for  cases  1-5  in  Figure  17  be 
explained  by  our  exclusion  of  diffusion  processes?  B^Jow  ~ 100  km 
eddy  diffusion  tends  to  dominate  transport  processes.^  For  an  eddy 
diffusion  coefficient,  D , approximately  independent  of  altitude  particles 
originally  at  a point  would  scatter  as  in  Brownian  motion; ^^nd,  over 
layers  of  thickness  d,  the  eddy  diffusion  time  is  given  by 

T = d^/2D^.  (9) 

“ 6 2 

For  d = 7.5  ^ (approximate  scale  height)  and  D = 10  cm  /s  we  find 
T = 2.8  X 10^  seconds.  This  computation  indicates  that  diffusion 
effects  are  unimportant  at  least  up  to  lO^s  and  that  the  build  up  of 
the  nitric  oxide  density  would  persist. 


gorj^the  nighttime  case,  Figure  14,  a quiescent  value  for  D is 
°s  and  we  see  that  none  of  the  curves  has  reached  this  value 
model  seconds. 


by 


j From  Equation  (8)  a detachment  process  that  is  important  under  all 

excitation  conditions  studied  is 

O'  + N2  N2O  + e R124 

This  reaction  had  been  assigned  the  upper  limit^*^  of  1 x 10  ^^cm^s’^. 

To  examine  the  effect  this  value  mi^^t  have  on  tlji|  ef|ec|ive  electron 
detachment  rate  we  assigned  a value  of  1 x 10  ^ ernes'  and  exercised 

the  model  for  cases  4 and  6 at  late  times  (10°-10^s).  The  results  are 
shown  in  Table  1.  For  case  6 there  is  less  than  20%  difference  in  the 
values  for  D for  the  two  values  of  *^229 • case  4 there  are  slightly 

^^See  for  example,  T.  Ogawa  and  T.  Shimazaki,  "Diurnal  Variations  of  Odd 
Nitrogen  and  Ionic  Densities  in  the  Mesosphere  and  Lower  Thermosphere: 
Simultaneous  Solution  of  Photochemical-Diffusive  Equations,"  J.  Geophys. 
Res.  3945-3960,  1975. 

12 

T.  Shimazaki  and  A.  R.  Laird,  "Model  Calculation  of  the  Diurnal 
Variation  in  Minor  Neutral  Constituents  in  the  Mesosphere  and  Lower 
Thermosphere  Including  Transport  Effects,"  J.  Geophys.  Res.  75, 
3221-3235,  1970. 

13 

F.  S.  Johnson  and  E.  M.  Wilkins,  "Thermal  Upper  Limit  on  Eddy  Diffusion 
in  the  Mesosphere  and  Lower  Thermosphere,"  J.  Geophys.  Res.  70,  1281- 
1284,  1970. 

14 

W.  Lindinger,  D.  C.  Albritton,  F.  C.  Fehsenfeld  and  E.  E.  Ferguson, 
"Reactions  of  O'  with  N2,  N2O,  SO2,  NH,,  CH.  and  C-H.  and  C,H,'  with 
0 from  300“ K to  Relative  Kinetic  Energies  of  ~ 2 eV,"  J.  ChemT  Phys. 

«,  3238-3242,  1975. 

^ F.  C.  Fehsenfeld,  private  communication,  1977. 
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TABLE  1.  COMPUTED  VALUES  OF  D (s‘^)  AS  A FUNCTION  OF  TIME  AND  VALUE 


OF  kj29  for  case 

6 AND  CASE 

4 

CASE 

10° 

10^ 

10^ 

loi 

10^ 

k 

'^129 

V time 

(cc/s) 

\(s) 

6 

10'^^ 

8.89(-l)** 

1.06(0) 

9.17(-1) 

8.60(-l) 

8.97(-l) 

6 

7.58(-l) 

8.95{-l) 

7.77(-l) 

7.50(-l) 

7.84(-l) 

4 

10'^^ 

8.14(+l) 

3.38(+l) 

4.98(0) 

9.77(-l) 

3.23(-2) 

4 

5.58(+l) 

2.27(+l) 

2.88(0) 

7.61(-1) 

2.92(-2) 

= 10^ 

-3  - 

ion  pairs  cm  s 

z 

o 

II 

o 

00 

-3 

cm 

*Qo 

= lo" 

-3  - 

ion  pairs  cm  s 

1 K1  mil  -3 

, Nq  = 10  cm 

**  -1 
Read  as  8.89  x 10  . 


) 
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differing  rates  of  decline  for  Jhe  two  values  of  kj-g  resulting  in  a 
maximum  of  53%  difference  at  10  seconds.  As  we  have  found  previously, 
with  the  smaller  value  for  k.-q*  reaction  129  was  no  longer  dominant,  but 
was  still  significant  (>  10%j  and  could  not  in  general  be  ignored. 

The  effective  recombination  coefficient  is  defined  by 

= (1  + + Xo.),  C9) 

where  X is  the  ratio  of  the  total  negative  ion  density  to  the  electron 
density  and  o.  is  the  effective  positive  ion-negative  ion  recombination 
coefficient.  In  the  steady  state,  a criterion  many  times  employed  in 
hybrid  models,^®  (solid  lines  in  Figures  18  and  19),  shoi^ld  equal  the 
source  term  divided  by  the  electron  density  squared,  Q/[e]  (dashed 
lines  in  Figures  18  and  19).  This  equality  is  approximately  reached 
(i^e.,  within  20%)  in  both  Figures  18  and  19  by  10”^s  for  curve  by 
10  s for  curves  2 and  3 and  by  lO'^s  for  curves  4,  5 and  6 (by  lO'^s  for 
curve  6 at  night) . These  curves  tend  to  be  grouped  according  to  the 
value  of  0^.  We  have  not  been  able  to  find  a simple  explanation  as  to 
why  this  is  so. 


The  governing  equation  for  the  approach  to  the  steady  state  is 

d[e]/dt  = Q(t)  - i|;{a^(Q),  X(Q),  a.(t)}[e]^  , (10) 

where  Q and  iji  are  given  by  equations  (1)  and  (9),  respectively.  The 
explicit  dependencies  are  written  out  to  emphasize  the  fact  that  li/  is 
an  implicit  function  of  time  through  Q.  From  equation  (10),  one  might 
expect  the  tendency  to  group  the  time  at  which  a steady  state  condition 
is  reached  with  Q , as  observed  in  Figures  18  and  19.  But,  the  non-linearity 
and  implicit  func?ional  dependencies  upon  time  in  equation  (10)  have 
prevented  any  simple  qualitative  understanding  of  the  observed  grouping 
with  0 . Several  attempts  to  solve  simpler  formulations  of  equation  (10) 
indic^e  that  these  simpler  solutions  would  not  reproduce  the  computed 
results  of  Figures  18  and  19.  In  other  words  to  attain  the  results  observed 
in  Figures  18  and  19,  we  have  found  that  equation  (10)  should  be  solved 
in  its  entirety,  as  is  accomplished  by  the  BENCHMARK-76  program. 


SUMMARY 

The  BENCHMARK-76  code  has  been  employed  to  obtain  species  densities 
as  a function  of  time  and  excitation  conditions  for  an  altitude  of  60 
km.  The  excitation  conditions  employed  were  as  follows.  The  prompt 

ionization,  N^,  was  assigned  the  values  10^\  loJ^2°^  cm”^,  and  the 

delayed  ionization  is  given  by  Q(t)  = Q^(l  + t)~  ’ where  t is  the  time 


See  for  example,  A.  P.  Mitra,  "D-Region  in  Disturbed  Condition, 
Including  Flares  and  Energetic  Particles,"  J.  Atmos,  and  Terr.  Phys. 
37,  895-913,  1975. 
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in  seconds  and  Q assigned  the  values  of  10^,  10®  or  10^^  ion-pairs 

cm"  s’  , subject  to  the  constraint  |q  | < |n  |.  Computations  were  made 
for  three  hour  intervals  around  noon  ^DAY)  and  midnight  (NIGHT) . 

Analysis  of  these  calculations  is  summarized  as  follows. 

First,  for  a given  delayed  ionization  parameter,  Q , it  is  found 
that  the  greater  the  initial  ionization  parameter,  N , ?he  sooner  the 
electron  density  relaxes  to  its  quiescent  value.  This  inversion  of 
one's  intuitive  feeling  is  traced  to  severe  neutral  chemistry  changes 
that  affect  the  negative  ion  chemistry. 

Second,  the  comparison  between  the  BENCHMARK-76  (BRL)  and  the  DCHEM 
(MRC)  calculations  of  the  electrgn  densities  at^bOj^km  for  both  ^ay  and 
ni^^t  conditions  and  for  0 = 10°  ion-pairs  cm  °s’  and  N^  = 10  or 
10^  cm"  are  good  to  exc^lent. 

Third,  the  calculated  relaxation  of  the  effective  electron-ion 
recombination  coefficient  to  quiescent  values  is  dominated  by  Q^,  i.e., 
the  larger  Q the  longer  the  relaxation  time.  The  return  to  quiescent 
values  for  tRe  higher  excitation  conditions  is  computed  to  take  longer  than 
10^  seconds. 

Fourth,  variations  in  the  effective  electron  attachment  frequency, 

A,  amounting  to  almost  a factor  of  two  at  night,  are  traced  in  detail 
and  found  to  result  from  ozone  enhancement  at  late  times. 

Fifth,  only  the  daytime  effective  electron  recombination  coeffi- 
cient for  the  mildest  excitation  studied  is  found  to  relax  to  a typical 
quiescent  value.  In  the  daytime,  all  other,  stronger  excitation  cases 
indicate  a relaxation  to  values  much  below  that  expected  for  quiescent 
conditions.  This  phenomenon  is  traced  to  the  greater  build  up  in  the 
nitric  oxide  density  with  the  larger  value  of  N^.  In  the  nighttime  ^ 
cases,  no  curves  approach  their  quiescent  values  even  after  10  or  10  s. 

We  have  found  that  the  reaction  O'  + N2  -►  N2O  + e is  non-negligi^^e,  ^ 
even  when  its  rate  coefficient  is  reduced  one  hundredfold  to  10  cm  /s. 

Finally,  the  time  |o  attain  a steady  state  condition  as  measured  by 
the  criterion  = Q/[e]  , is  found  to  be  a function  of  Q . The  smaller 
this  parameter,  the  longer  it  takes  to  reach  the  steady  state.  Attempts 
to  understand  this  behavior  by  simplifying  the  governing  non-linear 
differential  equation  for  the  electron  density  were  not  successful,  but, 
led  us  to  the  conclusion  that,  in  order  to  obtain  these  results,  this 
equation  ought  to  be  solved  in  its  entirety,  as  is  accomplished  by  the 
BENCHMARK- 76  program. 
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APPENDIX 


The  table  in  this  Appendix  displays  the  nominal  496  reactions  that 
link  the  64  species  within  the  BENCHMARK-76  code.  Reactions  proceed  from 
left  to  right  only.  Each  rate  coefficient,  k,  is  constructed  from  the 
reaction  parameters  A,  B and  C,  i.e. 

k = A(T/300)®  exp  (-  C/T)  , 

where  T is  the  temperature  in  degrees  Kelvin.  Except  for  the  first  60 
reactions,  involving  a photon  as  a reactant,  the  reference  for  the  value 
of  the  rate  coefficient  is  given  to  the  right  of  the  parameter  C column. 
Further  details  can  be  obtained  from  Reference  3. 
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